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The magnetic properties of the diluted magnetic semiconductors (DMS) (Ga, Mn)As and (Ga, 
Mn)N are investigated by means of an effective Heisenberg model, whose exchange parameters 
are obtained from first-principle calculations. The finite-temperature properties of the model are 
studied numerically using a method based upon the Tyablikov approximation. The method properly 
incorporates the effects of positional disorder present in DMS. The resulting Curie temperatures for 
(Ga, Mn)As are in excellent agreement with experimental data. Due to percolation effects and 
noncollinear magnetic structures at higher Mn concentrations, our calculations predict for (Ga, 
Mn)N very low Curie temperatures compared to mean-field estimates. 
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I. INTRODUCTION 



Ferromagnetic (III, Mn)-V diluted magnetic semicon- 
ductors (DMS) have attracted considerable attention 
among scientist during the past yeargi^. Their investi- 
gation has been driven by the idea of using their coupled 
electronic and magnetic degrees of freedom to construct 
electronic devices ranging from fast nonvolatile memo- 
ries to quantum computers 3 . To date, however, techni- 
cal applicability has been limited by the fact that most 
known DMS have Curie temperatures Tq below room 
temperatur o 2 i 4 i 5 i 6 i 7 . 

For the development of ferromagnetic DMS with higher 
Curie temperatures, it is important to understand theo- 
retically the magnetism in these materials and to develop 
theories which provide reliable qualitative and quantita- 
tive predictions. The magnetism in these materials is 
due to magnetic moments localized at magnetic impuri- 
ties, which interact with each other indirectly via holes 
in the valence and impurity band of the host semicon- 
ductor. Therefore, for the description, one often employs 
an effective Heisenberg model, whose exchange parame- 
ters are determined by the interaction between the local- 
ized moments and the holes^iiSiiiiiSiiiiiii 5 .. However, 
the magnetic impurities are mainly randomly distributed 
over the sites of the crystal lattice. This positional dis- 
order breaks the translational symmetry of the crystal 
and thus greatly complicates the theoretical description 
of the material. Studies based on the mean-field ap- 
proximation (MFA)^ or the random-phase approxima- 
tion (RPA) combined with the virtual-crystal approxi- 
mation (VCA)iS neglect effects of the positional disorder 
in DMS. Approaches based on percolation theor; 
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count for the randomness of the impurity positions, but 
require a simple functional dependence of the exchange 
parameters on the inter-spin distance and treat the mag- 
netism itself only on a mean-field level. Monte-Carlo sim- 
ulations (MC)^ 3 ^^^ 6 " seem to provide a better way to 
include the positional disorder, but these are numerically 



expensive and usually assume classical spins. However, 
a proper treatment of the positional disorder of the lo- 
calized moments and their quantum nature is needed to 
make reliable predictions about the magnetic properties 
of DMSiLi^. 

In a previously published article^, the exchange pa- 
rameters of an effective (classical) Heisenberg Hamil- 
tonian have been calculated from first-principles for 
Ga 1 _ a ,Mn 2 ,As and Ga 1 _ 2 ,Mn x N. There, however, these 
had only been used to calculate Curie temperatures 
within MFA. More recently, results of classical MC 
simulations on the basis of these exchange parameters 
have been presented^. Here, we employ a different 
approachi2i2& to investigate the properties of the effective 
Heisenberg Hamiltonian. This approach generalizes the 
Tyablikov approximation 21 to systems with positional 
disorder, which is treated numerically exactly. Further- 
more, the method assumes quantum spins. The quan- 
tum fluctuations of the spins are treated within random- 
phase approximation, which goes beyond MFA and the 
classical-spin approximation. It should be mentioned 
that a similar approach has been proposed in2£. 



II. MODEL 

Details of the electronic-structure calculation for 
Ga 1 _ 2 .Mn 2 ,As and Ga 1 _ 2 ,Mn ;E N and the extraction of the 
exchange parameters J(R) as function of the Mn-Mn dis- 
tance R can be found in RefiS. Here, these exchange 
parameters are used as input for a 'diluted' Heisenberg 
model 
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in which only a fraction of the lattice sites is occupied by 
a spin. Hence, i and j label the occupied lattice sites only, 
whose total number is N, and e i = (Sf,SV,Sf) /{hS) 
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is the normalized spin operator of the localized mag- 
netic moment at lattice site i with lattice vector R^ and 
= J(|R ; — R-yD- The magnitude S of the spins is ab- 
sorbed by the exchange parameters due to the particu- 
lar way in which these are calculated from the electronic 
structure. 

The finite-temperature properties of Hamiltonian 
are studied using a generalization of the Tyab- 
likov approximation to systems without translational 
symmetrjii2iSii. The generalization treats the positional 
disorder in the spin system numerically exactly except 
that a uniform magnetization is assumed. Further- 
more, the effects of low-energy quantum excitations, i.e. 
magnons, are included. Within this approximation, the 
local magnon spectral density is given byl&2&: 
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where the E are the eigenvalues of the Hamilton 
matrix H, which is defined by its matrix elements 



H i 3 = S ij S n =i J m ~ J ij- These eigenvalues also deter- 
mine the Curie temperature: 
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To evaluate this expression for a given set of E r 's, the 
value of S has to be fixed. For Mn ions in Ga 1 _ 2 ,Mn 2 ,As 
and Ga 1 _ x Mn a .N, S = 5/2 should be appropriate 2 . How- 



ever, this choice is not consistent with the calculation of 
the exchange parameters from the electronic structure, 
where classical spins are assumed. Therefore, we will use 
Eq. © in the limit S — > 00, which yields T c values a 
factor 5/7 less than for S = 5/2. 

Due to the positional disorder of the spins present in 
DMS, the eigenvalues cannot be computed by Fourier 
transformation of H. However, the eigenvalues may be 
obtained by the numerical diagonalization of the Hamil- 
ton matrix for a finite system. In our calculations, we 
used systems of ~ 10 000 spins, which were randomly 
distributed over the lattice sites of a cubic section of an 
fee lattice with periodic boundary conditions. For each 
concentration x of Mn ions, we averaged the spectral den- 
sities over eight random configurations. 



III. RESULTS 



In Fig. the Mn-Mn exchange interactions J(R) in 
Ga, Mn As and in Ga x _ Mn N are shown as functions 
of the Mn-Mn distance R for several concentration x. In 
Ga x Mn As, the falloff of the interaction with R is com- 
parably slow. In Ga 1 _ 2 ,Mn 2 .N, the interaction between 
nearest neighbors is much larger than in Ga 1 _ a ,Mn a; As, 
but Mn moments further apart are only very weakly cou- 
pled. 



FIG. 1: Exchange interactions J(R) between Mn ions of dis- 
tance R in (a) Gai-^Mn^As and (b) Gai-^Mn^N for various 
concentrations x (fror 9 '" ' 1 



Figure [21 shows the resulting magnon spectral densi- 
ties. For Ga 1 _ x Mn 2 ,As, the spectrum is smooth and 
continuous. For Ga 1 _ x Mn x N, one can recognize rem- 
nants of peaks typical for nearest-neighbor interaction at 
low concentrations, which are broadened by small long- 
ranged interactions. Compared to Ga 1 _ x Mn x As, there is 
a large spectral density at low energies for Ga^ Mn N. 
For concentrations x > 0.08, antiferromagnetic interac- 
tions come into play and negative magnon energies ap- 
pear indicating a ground state which is not a saturated 
ferromagneliSil. 

The Curie temperatures calculated using Eq. are 
shown in Fig. For Ga 1 _ 3 ,Mn a .As, the calculated val- 
ues agree remarkably well with the experimental values 
of optimally annealed samples^24. Furthermore, the 
calculated curve suggests that slightly higher T c 's might 
be achieved by further increasing the Mn content x, but 
values above 300K seem rather unlikely. 

Since experimental values for T c in Ga 1 _ x Mn x N are 
quite controversial (reported values range from OK to 
940K 4 - 26 - 27 *22i2iL2ii), we refrain from a comparison here. 
However, the Curie temperatures we calculated are quite 
low compared to earlier mean-field estimates (e.g. in 
Refi£). These low T c values despite the high values of 
the nearest-neighbor exchange may be explained as fol- 
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FIG. 2: Local magnon spectral density Sa(E) for (a) 
Gai-jjMriajAs and (b) Gai-^Mn^N for various concentrations 
x of Mn 



FIG. 3: Calculated Curie temperature Tc of (a) 
Gai-jMnjAs (compared with experimental values of an- 
nealed samples"*"*—*"-) and (b) Gai-^Mn^N for various con- 
centrations x of Mn ions 



lows: For concentrations well below the nearest-neighbor 
percolation threshold c P « 0.2*2i, even a large nearest- 
neighbor exchange does not contribute substantially to 
the stability of the magnetic phase. Since the exchange 
parameters for larger inter-spin distances are very small 
in Ga-, Mn N, ferromagnetic order can only be estab- 
lished at very low temperatures. Note that the drop of 
T c for x > 0.08 may be due to the used approximation: 
As indicated by the magnon spectra seen in Fig. [2 the 
system's ground state is different from a saturated ferro- 
magnet, but a such uniform magnetic state is assumed in 
the approximation. 

Figure ^ presents a comparison of the the Curie tem- 
peratures calculated using different approximations for 
the effective Heisenberg model. The T c values ob- 
tained by MC simulations are slightly higher than the 
ones calculated by the presented approach, whereas 
both MFA and VCA-RPA yield much higher T c 's. For 
Ga 1 _ 2 ,Mn :z , As, the difference is about a factor two to 
eight. For Ga, Mn N, the difference is even much 
larger. This is due to the fact that the MFA and VCA- 
RPA do not take into account percolation effects: Large 
nearest-neighbor interactions yield large Curie temper- 
atures even for concentrations well below the nearest- 
neighbor percolation threshold. However, for such con- 
centrations, the nearest-neighbor interaction strength 



should not play an important role for the ferromagnetic 
stability, which can be easily seen by considering the case 
of nearest-neighbor interaction onljiSS. 

IV. SUMMARY 

In this paper, we presented a method for calculat- 
ing the magnetic properties of ferromagnetic DMS. The 
method applies a Tyablikov-like approximation for sys- 
tems with positional disorder to an effective Heisenberg 
Hamiltonian, whose exchange parameters where obtained 
by first-principle calculations. Unlike in MFA or VCA- 
RPA, no approximations with respect to the positional 
disorder are made apart from the simplification of a uni- 
form magnetization. As the main advantage over clas- 
sical MC simulations, the presented treatment of the ef- 
fective Heisenberg model admits quantum spins and thus 
may open up a way towards a fully quantum-mechanical 
treatment of magnetism in DMS. Furthermore, the nu- 
merical effort is fairly low compared to MC simulations. 

Our calculations of T c for Ga 1 _ x Mn As show excellent 
agreement with experimental data. For Ga, Mn N, 
we obtained very low Curie Temperatures despite high 
effective nearest-neighbor exchange parameters, which 
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FIG. 4: Comparison of the Curie temperatures Tc of 
Gai-jMnjAs (diamonds) and Gai-^Mn^N (squares) ob- 
tained by the presented approach (solid line, filled symbols), 
VCA-RPA (dashed line, filled symbols), MFA (dotted line, 
filled symbols) and MC (dash-dotted line, open symbols, 
taken from—). 



shows the importance of percolation effects. Moreover, 
for both Ga, Mn As and Ga, Mn N, the T n values 

_L — X X _L — X X ' O 



we found are much lower than MFA and VCA-RPA val- 
ues. These results support recent findings obtained by 
using MC simulations in combination with first-principle 
methods 1 ^. 

The presented model should be improved by using a 
self-consistent method describing the electronic degrees 
of freedom at finite temperature (such as, e.g., in— iS). In 
order to obtain a fully quantum mechanical theory, quan- 
tum spins should be used instead of classical spins in in 
the calculation of the effective exchange parameters from 
the electronic structure. This will also remove the ambi- 
guity in the choice of S. Furthermore, the treatment of 
the effective Heisenberg model may be extended to allow 
for a site-dependent (Sf). In addition, the model might 
be improved in order to handle systems with a ground 
state deviating from a saturated ferromagnet. Further- 
more, clustering and other forms of short-range chemical 
ordering may also be included into the model in order to 
investigate their effects on the magnetic stability. Finally, 
the method should be applied to other DMS. 
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